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Abstract 

We study the regime of anticipated synchronization in unidirectionally coupled 
chaotic maps such that the slave map has its own output reinjected after a cer- 
tain delay. For a class of simple maps, we give analytic conditions for the stability 
of the synchronized solution, and present results of numerical simulations of cou- 
pled ID Bernoulli-like maps and 2D Baker maps, that agree well with the analytic 
predictions. 
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The synchronization of chaotic systems is a subject that has attracted a lot 
of attention in the past years. Since the pioneering works [1] several different 
regimes of synchronization have been found: complete synchronization, phase 
synchronization [2], lag synchronization [3], generalized synchronization [4,5], 
synchronization by common noise force [6,7], among others. 

Anticipated synchronization is a recently discovered synchronization regime 
that occurs in unidirectionally coupled systems [8,9]. In this regime coun- 
terintuitive phenomena occur, since the slave system anticipates the chaotic 
evolution of the master system, despite the fact that chaotic behavior implies 
long-term unpredictability. 
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In the case of coupled time-delayed differential equations the anticipation time 
is related to the difference between the intrinsic delay time of the systems and 
the delay time of the coupling [10]. In coupled ordinary differential equations 
(as the Lorenz and Rossler systems) the anticipation time must be small to 
have a stable synchronization manifold [8] . However, by using a chain of slave 
systems, anticipation times that are multiples of the coupling delay time and 
that exceed characteristic time scales of the chaotic dynamics can be obtained 
[9,11]. Many additional numerical [12,13] and experimental [14,15] studies of 
anticipated synchronization have been performed. 

Recently, analytic conditions for the synchronization of coupled maps with 
delays were given by Masoller and Zanette [16]. In the maps considered in [16] 
the chaotic behavior is induced by the delay term in the map (in other words, 
without delay the master map is not chaotic). In Ref. [16] the master map 
(x n ) and the slave map (y n ) are of the form 

x n+1 = bx n + f(x n - N ), (1) 

Vn+i = by n + (1 - r})f(y n _ N ) + vf( x n-M)- (2) 

with \b\ < 1 and / a nonlinear function. The synchronized solution is given by 
Vn — Xn-M+N- By studying the evolution of a small perturbation of the solu- 
tion it was shown that the synchronization is stable for r] — 1. The existence 
of a threshold value of the coupling, r\ c < 1, above which the synchronized 
solution is stable, was also shown. 

When the chaotic behavior of the master map is not induced by a delay, 
to the best of our knowledge no analytical conditions for anticipation have 
been reported. In this letter we study a system composed by a chaotic map 
(master) unidirectionally coupled to a second chaotic map (slave) which has 
its own signal reinjected after a certain delay. We consider a master map, x n , 
and a slave map, y n of the form 

x n+ i = f(x n ), (3) 

Vn+l = f{Vn) + l(Xn-N ~ Vu-m), (4) 

where / is a nonlinear function. We show that the delayed auto-injection in 
the slave map leads to anticipation in the synchronization, and present ana- 
lytic conditions for the stability of anticipated synchronization. An analytical 
treatment is possible because we consider simple chaotic maps. We exemplify 
the results with numerical simulations of coupled ID Bernoulli-like maps and 
2D Baker maps. 
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The synchronization manifold of (3), (4) is given by 



y n = x n - N+M , (5) 



and thus the slave variable is lagged by N — M steps behind the value of the 
master variable (if N — M < 0, the slave map anticipates the dynamics of the 
master map). To study the stability of the synchronized solution we consider 
a perturbation of the form 

y n = x n - N+M + S n . (6) 



In the linear regime, the perturbation obeys the following map 

(7) 

where /' = df/dx. Making the change of variables 8 n = (z n ) n gives 

(z n+1 ) n+1 = f(x n . N+M )(z n ) n - 7(^-M) n_Af . (8) 



Notice that the value of N is irrelevant in the long-term behavior, thus the 
stability of the synchronized solution is the same in the anticipated and in the 
retarded regimes (similar results were found in [16]). Therefore, and without 
loss of generality, in the following we consider N = 0. 

A sufficient condition for the stability of the synchronized solution will be that, 
for n > no, where no is some number of transient steps, all the solutions z l n of 
(8) (where % labels the different solutions) satisfy \z\\ < 1. When /' depends 
on x n , an analytical treatment of the stability of the synchronized solution is 
in general not possible. However, there are particular cases in which the study 
of Eq. (8) gives insight into the parameter region where synchronization is 
stable. As a first example we consider the ID Bernoulli-like map 

f(x n ) = ax n mod 1 . (9) 



The map is chaotic for a > 1. /' = a, and thus the solutions of (7) are linear 
combinations of functions of the form S n = z n , with constant z. Thus (8) reads 

z M+1 = az M - 1 . (10) 



Stability of the synchronized solution is obtained if all the solutions of (9) 
satisfy \z\ < 1. 
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For M = 1 the roots of Eq. (10) are the solutions of a simple quadratic 
equation, so that it is simple to check that they have \z\ < 1 when 7 G (a— 1, 1). 
For arbitrary M a necessary (but not sufficient) condition for the stability of 
the synchronized solution is 7 G [a — l,a + 1]. To show that, first note that 
for 7 = 0, Eq. (10) has one root at z = a > 1, and the other M roots are 
at z = 0. A simple analysis of perturbations shows that a small 7 breaks the 
degeneracy of the roots at the origin, which then move radially outwards in 
the complex plane, whereas (if 7 > 0) the root that was located at z = a 
diminishes its value, moving towards the unit circle. Thus, stability of the 
synchronized solution will be obtained by increasing 7 if this last root enters 
into the unit circle in the complex plane before some of the other M roots 
leave it. Finally, for 7 large enough, all the roots are outside the unit circle. 
Eq. (10) can be rewritten as 

z M = 1 —. (11) 



Thus, the value of \z\ is bounded between 



< < \ z \m < —L—. (12) 



a + \z 



At the limits of the stability region, some root z would satisfy \z\ = 1, so that 

< 1 < (13) 
a + 1 ~ ~ a- 1 v ; 



which leads to a — 1 < 7 < a + 1. Therefore, only within this interval a root 
might cross the unit circle, leading to stability changes of the synchronized 
solution. In consequence the range of 7 leading to stable synchronization is 
inside 7 G [a — 1, a + 1]. Further insight can be gained from the study of the 
roots at the boundaries of this interval. First we consider the case 7 = a — 1. 
For this value of the coupling Eq. (10) becomes 

2"+! = az n -a + l. (14) 



Clearly, z = 1 is a solution for all a and M. A small perturbation of the value 
of 7, 7 = a — 1 + ^7 leads to a modification of the value of z, z — 1 + Sz. To 
first order, 5z and £7 are related by 

5z = 5 4 r. (15) 

l-M(a-l) v ; 

The denominator is positive if M < l/(a — 1) and negative otherwise. Thus, 
if M < l/(a — 1), by increasing 7, (#7 > 0), Sz < and a real root enters 
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into the unit circle. Since this value of 7 is the smallest one for which crossing 
the unit circle becomes possible, all the roots satisfy now \z\ < 1 and thus the 
synchronized solution becomes stable. On the other hand, if M > l/(a — 1), 
by increasing 7, 5z > 0, and a real root leaves the unit circle. In this case the 
synchronized solution becomes more unstable, in the sense that the rate of 
escape given by the largest \z\ increases by increasing 7. Thus, M < l/(a — 1) 
gives a limit of the number of steps M for which the synchronized solution 
can become stable. There is a relation between M and the degree of chaos, 
associated to the Lyapunov exponent log a, of the master map: the largest the 
value of a, the lower the stable anticipation times. If a > 2, stable anticipated 
synchronization becomes impossible. 

Next, we consider the other boundary of the synchronization region, 7 = a + 1. 
In this case Eq. (10) becomes 

z** 1 = az n -a-l. (16) 

If M is even, z = — 1 is a solution. Considering a perturbation of the form 
7 = a + 1 + ^7, in the same way as before it can be shown that z = — 1 + 5z 
with 5z = —Sj/ll + M(a + 1)]. Since the denominator is always positive, if 
7 grows, z decreases, so that the root close to z — — 1 leaves the unit circle. 
Thus, we simply confirm that a necessary condition for the stability of the 
synchronized solution is 7 < a + 1. 

Next we show results of numerical simulations that confirm these analytic 
arguments. Figure 1 shows simulations of Bernoulli-like maps, in which the 
slave anticipates the master in five steps. During the first 5000 steps the maps 
evolve independently. Then we set the value of the coupling to 7 = 0.15 
and after a very short transient the anticipation of the slave to the master 
is evident. Figure 2 shows, for the same value of the parameters a and M, 
and increasing coupling 7, how the roots of Eq. (10) move in the complex 
plane. For low coupling Eq. (10) has one real root z\ ~ a > 1, one real root 
Z2 < 1, and four complex conjugate roots with modulus less than 1 [Fig. 2 
(a)]. As the coupling increases, z\ decreases while the other roots increase 
their modulus, approaching the unit circle. For 7 > 0.1 = a — 1 all roots of 
Eq. (8) have modulus less than 1 [Fig. 2 (b,c)] and the synchronized solution 
is stable (Fig. 2(c) corresponds to Fig. 1). For even larger coupling, pairs of 
complex-conjugate roots cross the unit circle [Fig. 2 (d)], and synchronization 
is unstable again. 

For values of the coupling such that all roots of Eq. (10) are inside the 
unit circle, the distance between the two trajectories decreases exponentially, 
\x h+ m — Un\ ~ \xm — Uo I exp(— n/r), with the transient time to synchroniza- 
tion given by the inverse of the logarithm of the modulus of the largest root, 
r = — 1/ In l^i I . Fig. 3 shows the transient evolution of \x n+ M — y n \, for the 
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Fig. 1. Time series of Bernoulli-like maps (x n solid line, y n dashed line) for a — 1.1, 
7 = 0.15, M = 5. The coupling is set on at n = 5000, leading to anticipated 
synchronization. 

parameters of Fig. 1. We observe a damped oscillatory behavior. While the 
damping time is r, the frequency is associated to the phase of z\. 

As a second example we consider the 2D Baker's map that transform the unit 
square into two non overlapping rectangles: The master map, and 
the slave map, are 



~m _ f ( m m\ 

C + i = /,(CC) 



(17) 
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Fig. 2. Roots of Eq. (10) for a = 1.1, M = 5 and (a) 7 = 0.005, (b) 7 = 0.11, (c) 
7 = 0.15, and (d) 7 = 0.8. 



<+i = f x (x s n , y s n ) + t(C - K-m) 

Vn+1 = Vn) + TfiVn ~ Vu-m)- 



(18) 



where 



fx 



fy = < 



ax n if x n < 1/a, 

a(x n — 1/a) if x n > 1/a, 

by n tix n <l/a, 
by n + (1 - b) if x n > 1/a, 



(19) 



(20) 



a > 1 and b < 1 are the expansion and contraction rates, respectively. It 
is easy to see that for the stability of the synchronized solution we obtain a 
pair of equations of the same form as Eq. (10), where now a is equal to the 
expansion and to the contraction rate, respectively. Both equations must have 
roots with modulus less than 1, for the synchronized solution to be stable. 
In Fig. 4 we present results of numerical simulations that show anticipated 
synchronization by one step. The coupling is set on at n = 30, and after a 
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Fig. 3. Transient decay of |x„+m — Un\ for a = 1.1, M = 5 and 7 = 0.15. The solid 
line indicates the value of |x n +M — Un\ and the dashed line indicates the value of 
\xm — 2/0 1 exp(— n/r) with r = — 1/ln \z\ \ ~ 15.55 as explained in the text. 



transient the map solution approaches the anticipated synchronization state. 
We find synchronized solutions for parameter values such that the equations 
analogous to (10) have roots with modulus less than 1. The duration of the 
transient, again, is related to the modulus of the largest root. 

In summary, we have studied the regime of anticipated synchronization in uni- 
directionally coupled chaotic maps. In a general case is not possible to give 
analytic conditions for the parameter region where anticipated synchroniza- 
tion occurs, but we have presented two class of maps, ID Bernoulli-like maps 
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Fig. 4. Time series of the 2D Baker map (x™ solid line, x s n dashed line) for a = 1.333, 
b = 0.777, 7 = 0.7, and M = 1. The coupling is set on at n = 30, leading to 
anticipated synchronization. 



and 2D Baker maps, in which an analytic treatment of the stability of the 
synchronized solution is possible. The results of numerical simulations are in 
good agreement with the analytic predictions. 
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